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ABSTRACT 

The integrated Spectral Energy Distributions (SED) of the Large Magellanic Cloud (LMC) and Small Magellanic Cloud (SMC) appear significantly 
flatter than expected from dust models based on their far-infrared and radio emission. The origin of this millimetre excess is still unexplained, 
and is here investigated using the Planck data. The integrated SED of the two galaxies before subtraction of the foreground (Milky Way) and 
background (CMB fluctuations) emission are in good agreement with previous determinations, confirming the presence of the millimetre excess. 
The background CMB contribution is subtracted using an Internal Linear Combination (ILC) method performed locally around the galaxies. The 
foreground emission from the Milky Way is subtracted as a Galactic H i template and the dust emissivity is derived in a region surrounding the two 
galaxies and dominated by Milky Way emission. After subtraction, the remaining emission of both galaxies correlates closely with the atomic and 
molecular gas emission of the LMC and SMC. The millimetre excess in the LMC can be explained by CMB fluctuations, but a significant excess 
is stiU present in the SMC SED. 

The Planck and IRAS-IRIS data at 100 |^m are combined to produce thermal dust temperature and optical depth maps of the two galaxies. The 
LMC temperature map shows the presence of a warm inner arm already found with the Spitzer data, but also shows the existence of a previously 
unidentified cold outer arm. Several cold regions are found along this arm, some of which are associated with known molecular clouds. The dust 
optical depth maps are used to constrain the thermal dust emissivity power law index (p). The average spectral index is found to be consistent with 
P =1.5 and jS =1.2 below 500 |xm for the LMC and SMC respectively, significantly flatter than the values observed in the Milky Way. Furthermore, 
there is evidence in the SMC for a further flattening of the SED in the sub-mm, unlike for the LMC where the SED remains consistent with =1.5. 
The spatial distribution of the millimetre dust excess in the SMC follows the gas and thermal dust distribution. 

Difl'erent models are explored in order to fit the dust emission in the SMC. It is concluded that the millimetre excess is unlikely to be caused by 
very cold dust emission and that it could be due to a combination of spinning dust emission and thermal dust emission by more amorphous dust 
grains than those present in our Galaxy. 

Key words. ISM: general, dust, extinction, clouds - Galaxies: ISM - Infrared: ISM - Submillimeter: ISM 



1. Introduction 

Star formation and the exchange and evolution of materials be- 

* Corresponding author; email: Jean-Philippe.Bernard@cesr.fr. tween the stars and the interstallar medium (ISM) are continuous 
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processes that drive the evolution of galaxies. As stars evolve, 
die and renew the life cycle of dust and gas, the amount of dust 
and its distribution in a galaxy has important consequences for 
its subsequent star formation. Thus knowing the dust content 
throughout cosmic history would provide clues to the star for- 
mation history of the universe as the metallicities evolve. One 
of the puzzling results that has emerged from the studies of dust 
SEDs in the early universe, is finding very high dust masses in 
high redshift galaxies. For example, recent Herschel observa- 
tions of submillimetre galaxies (SMGs) find excessively high 
dust masses, and high dust-to-gas mass ratios (D/G) (Santini 
et al. 2010). These findings contradict the low metallicities mea- 
sured in the gas. Large dust masses have also been measured for 
a range of low metallicity local universe dwarf galaxies (Dumke 
et al. 2003; GalHano et al. 2003, 2005; Bendo et al. 2006; 
Galametz et al. 2009; Grossi et al. 2010; O'Halloran et al. 2010). 
Our understanding of how dust masses are estimated and what 
observational contraints are necessary is called into question by 
the discovery that the dust masses measured in these galaxies ap- 
pear significantly higher than expected from their metal content. 

Excessively large dust masses only seem to be found in low 
metallicity galaxies to date. The evidence is found in the be- 
havior of the submillimetre (submm) emission, beyond about 
400 - 500 |Ltm. From a study of a broad range of metal-rich and 
metal-poor galaxies observed with wide wavelength coverage 
that included submm observations beyond 500 jam, a submm ex- 
cess, beyond the normal dust SED models, was found only for 
low metallicity systems (Galametz et al. 2010). Such submm ex- 
cess is particularly evident in the Magellanic Clouds (Israel et al. 
2010; Bot et al. 2010b). The origin of such excess is still uncer- 
tain but several suggestions have been put forward: 

- The submm excess has been modeled as a cold dust compo- 
nent with a submm emissivity index (fi) of yS=l, which sug- 
gests a low dust temperature of ~10K (Galliano et al. 2003, 
2005; Galametz et al. 2009, 2010). There are few observa- 
tional constraints in the submm wavelength window, hence 
this cold component is a rather ad hoc solution. This descrip- 
tion often results in discrepantly large D/G ratios being found 
for the low metallicity galaxies if the estimates of the total 
gas reservoir are well known. 

- Meny et al. (2007) show that emission by amorphous grains 
is expected to produce excess emission in the submm. As a 
result, the spectral shape of emission by amorphous grains 
cannot be reproduced by a single emissivity index over the 
whole Far-InfraRed (FIR) to millimetre (mm) range, and is 
expected to flatten at long wavelengths. The intensity of the 
excess is also predicted to depend strongly on the tempera- 
ture of the dust grains. Modifying the dust optical properties 
to incorporate the efl'ects of the disordered structure of the 
amorphous grains, showed that the efl'ective submm emis- 
sivity index decreases (thus flattening the submm Rayleigh- 
Jeans slope) as the dust temperature increases. Similarly, 
an inverse correlation between the dust temperature and the 
spectral index has been observed in data from the FIR to 
the submm (e.g. Dupac et al. 2003; Paradis et al. 2010). 
However, these authors as well as Shetty et al. (2009a) advise 
caution in the interpretation of the observed inverse T-p rela- 
tionship, which is afl'ected by the natural correlation between 
these two parameters, and requires careful treatment of ob- 
servational uncertainties. Shetty et al. (2009b) also studied 
the efl'ect on the observed T-p relation of temperature vari- 
ations along the line-of- sight. Paradis et al. (2009) using the 
DIRBE , Archeops and WMAP data showed that the FIR- 



submm SEDs of galactic molecular clouds and their neutral 
surroundings indeed show a flattening beyond /I ^ 500 |Ltm. 

- Spinning dust emitted in the ionised gas of galaxies has been 
suggested as the explanation for the radio emission often 
seen in galaxies (Ferrara & Dettmar 1994). This idea was 
further explored by Anderson & Watson (1993) and Draine 
& Lazarian (1998b,a) who characterized the 10 to 100 GHz 
anomalous foreground emission component as spinning dust 
from very small dust particles such as Polycyclic Aromatic 
Hydrocarbons (PAHs). Recent improvements to the model 
suggest the peak frequency, which could occur in the submm, 
depends on many parameters, including the radiation field 
intensity, the dust size distribution, the dipole moment dis- 
tribution, the physical parameters of the gas phases, etc. 
(Hoang et al. 2010; Silsbee et al. 2010; Ysard et al. 2010; 
Ysard & Verstraete 2010a). Spinning dust has been a pre- 
ferred explanation for the submm excess observed on the 
global scale in the Magellanic Clouds (Israel et al. 2010; Bot 
et al. 2010b). 

- ? have suggested hydrogenated amorphous carbon as the 
most likely carbonaceous grain material instead of graphite. 
The advantage of using amorphous carbon, particularly in 
the cases where a submm excess is found, is that amorphous 
carbons are more emissive and result in a flatter submm 
slope and less dust mass. This has been used with new 
Herschel observations to model some low metallicity galax- 
ies (Galametz et al. 2010; O'Halloran et al. 2008; Meixner 
et al. 2010). 

- The enhancement of hot, very small, stochastically heated 
grains with a low dust emissivity, has also been suggested 
to characterize the submm emission when there is a submm 
excess present (Lisenfeld et al. 2001; Zhu et al. 2009). 

Some of the possible causes of the submm excess could be 
ruled out or constrained if the gas mass estimates were better 
established. However, estimates of the total gas reservoir in low 
metallicity environments have been uncertain, particularly the 
molecular gas component. Observations of CO in low metallic- 
ity galaxies have been a challenge, and the dearth of detected 
CO has often been interpreted as meaning there is very little 
mass in molecular gas present in galaxies which are otherwise 
actively forming massive stars (Leroy et al. 2009). However, 
relying on CO alone to trace the H2 mass could potentially 
miss a large reservoir of molecular gas residing outside the CO- 
emitting region, as indicated by the FIR fine structure observa- 
tions in the Magellanic Clouds and other low metallicity galaxies 
(e.g. Poglitsch et al. 1995; Madden et al. 1997; Madden 2000). 
The presence of a hidden gas mass in low metallicity galaxies 
is also suggested by dust observations of the Magellanic Clouds 
(e.g. Bernard et al. 2008; Dobashi et al. 2008; Leroy et al. 2007; 
Roman-Duval et al. 2010; Bot et al. 2010a). 

Most of these studies which find a submm excess in low- 
metallicity galaxies are derived from models on global galaxy 
scales and use data with a wavelength coverage that is limited to 
X < 500 |Ltm with Herschel or to /I < 870 jam for ground-based 
observations. The Galametz et al. (2010) study demonstrated 
that dust masses can diff'er by an order of magnitude depending 
on the availability of submm constraints on the SED modelling. 
These require wider wavelength coverage of the submm and mil- 
limetre (mm) region and better spatial information, to map out 
the submm excess and determine the local physical conditions. 

The closest low-metallicity galaxies are our neighboring 
Large Magellanic Cloud (LMC) and Small Magellanic Cloud 
(SMC), which have metallicities of 50% and 20% solar metal- 
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licity (Zq) respectively. They are ideal laboratories for study- 
ing the ISM and star formation in different low metallicity en- 
vironments. Recent studies with Spitzer (SAGE: Meixner et al. 
2006; Bolatto et al. 2007) and Herschel (HERITAGE: Meixner 
et al. 2010) have mapped out the the temperature and dust mass 
distribution (Bernard et al. 2008; Gordon et al. 2009, 2010; 
Leroy et al. 2007). Recent Herschel studies of the LMC already 
point to a possible submm band excess and an excessive dust 
mass, which may be suggesting the presence of molecular gas 
not probed by CO. This could be explained by the presence of 
amorphous carbon grains (Roman-Duval et al. 2010; Meixner 
et al. 2010; Gordon et al. 2010). It was difficult to be conclusive 
about the cause of the excess seen in the Herschel observation 
because of the lack of longer wavelength coverage. Global ex- 
cess mm and submm emission in the LMC and SMC was in- 
disputably shown by Bot et al. (2010b) and Israel et al. (2010) 
using broader wavelength coverage that included the submm to 
cm observations of the Top-Hat balloon experiment, WMAP and 
COBE and their results favour spinning dust as the explanation 
for the submm and mm excess, but require additional spatial cov- 
erage to be firmly conclusive. Planck ^ observations have both 
wide wavelength coverage and the spatial resolution to locate 
the submm emission within the LMC and SMC, which will help 
determine the origin of the submm excess. 

2. Observations 

2.1. Planck data 

The Planck first mission results are presented in The Planck 
Collaboration (2011a) and the in-flight performances of the 
two focal plane instruments HFI (High Frequency Instrument) 
and LFI (Low Frequency Instrument) are given in The Planck 
Collaboration (201 le) and The Planck Collaboration (201 Id) re- 
spectively. The data processing and calibration of the HFI and 
LFI data used here is described in The Planck Collaboration 
(201 Ig) and The Planck Collaboration (201 Ih) respectively. 

Figure 1 shows the total intensity maps observed toward the 
LMC and SMC at a few HFI and LFI frequencies. Both galax- 
ies are well detected at high frequencies. Around 100 GHz, their 
emission can barely be distinguished from CMB fluctuations. At 
lower frequencies, the contrast between the galaxies' emission 
and the CMB fluctuations becomes larger again. Note that the 
apparent variation of the noise level at a constant right-ascension 
value across the LMC is real and is due to the LMC being posi- 
tioned at the edge of the Planck deep field. 

The Planck DR2 data have had the CMB fluctuations re- 
moved in a way which is inappropriate for detailed examination 
of some foreground sources like the LMC and SMC. We there- 
fore use the original data before CMB subtraction (version DX4 
of the Planck-HFI and Planck-LFI data) and perform our own 
subtraction of the CMB fluctuations, based on a local ILC (see 
Sec. 2.3.2). 

We use the internal variance (cTjj) provided with the Planck 
data, which represents the white noise on the intensity. The LMC 
and SMC have been observed many times, as the Planck scan- 
ning strategy (The Planck Collaboration 2011a) covers the re- 
gion close to the ecliptic pole repeatedly. In particular, the LMC 

^ Planck (http://www.esa.int/P/a;ic^) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries: 
France and Italy) with contributions from NASA (USA), and telescope 
reflectors provided in a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 



Table 1. Characteristics of the FIR/submm data used in this 
study. 
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* From Miville-Deschenes & Lagache (2005) 

^ Icr average value in one beam scaled from The Planck Collaboration 
(201 Ig). We actually use the internal variance maps for cfh 
^ From Bennett et al. (2003) 



is located at the boundary of the Planck deep field. Its eastern 
half has a much higher signal to noise ratio than its western half. 
We assume the absolute uncertainties from calibration given in 
The Planck Collaboration (201 Ig) and The Planck Collaboration 
(201 Ih) for HFI and LFI respectively and summarized in 
Table 1. 

2.2. Ancillary data 

2.2.1. H I LMC/SMC data 

To trace the atomic gas in the Magellanic Clouds, we used H i 
maps in the 21cm line, obtained by combining data from the 
Australia Telescope Compact Array (ATCA ) and the Parkes 
single dish telescope. For the LMC, this data was obtained by 
Kim et al. (2003) and Staveley-Smith et al. (2003) and covers 
11.1° X 12.4° on the sky. The spatial resolution is V correspond- 
ing to a physical resolution of about 14.5 pc at the distance of 
the LMC. For the SMC, the data was obtained by Staveley-Smith 
et al. (1997) and Stanimirovic et al. (1999). It covers a 4.5° x4.5° 
region. Hi observations of the SMC tail (7° x 6°) were obtained 
by Muller et al. (2003) and combined to the one in the direction 
of the SMC. The spatial resolution is 98'', corresponding to 30pc 
at the distance of the SMC. 

2.2.2. H I Galactic data 

The LMC and SMC are located at galactic latitudes bii=-34° 
and bn=-44° respectively. They can therefore suff'er from sig- 
nificant contamination by Galactic foreground emission, which 
has to be removed from the IR data. In order to account for the 
Galactic foreground emission, we used Galactic H i column den- 
sity maps. 

For the LMC, the Hi foreground map was constructed by 
Staveley-Smith et al. (2003) by integrating the Parkes H i data 
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Fig.l. Planck total Intensity data for the LMC (left) and SMC (right) at 857 (top), 100 (middle) and 28.5 GHz(bottom) at full 
resolution. The top panels are shown in log scale. The circle in the top panels shows the region used to extract average SEDs. 



in the velocity range from -64 < Vhei < 100 kms"\ which ex- 
cludes all LMC and SMC associated gas (v > 100 km s"^) but 
includes essentially all Galactic emission. The spatial resolution 
is 14'. For the SMC, a map of combined AJCA and Parkes data 
was build by integrating the Galactic velocities by E. Muller (pri- 
vate communication). The resolution is 98''. 



These maps show that the Galactic foreground across the 
LMC is as strong as A^h = 1.3 x 10^^ Hcm"^, with significant 
variation across the LMC, in particular a wide filamentary struc- 
ture oriented southwest to northeast. The Galactic foreground 
across the SMC is weaker with values A^h - 3.2 x 10^^ Hcm"^ 
but the associated FIR-submm emission is still non-negligible 
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since the SMC emission is also weaker than that of the LMC, 
due to its lower dust and gas content. 

These foreground maps are used to subtract the foreground 
IR emission from the IR maps, using the emissivity SED de- 
scribed in Sect. 2.3.3. 

2.2.3. CO data 

The ^^CO(/=1^0) molecular data used in this work was ob- 
tained using the NANTEN telescope, a 4-m radio telescope of 
Nagoya University at Las Campanas Observatory, Chile (see 
Fukui et al. 2008). The observed region covers about 30 square 
degrees where CO clouds were detected in the NANTEN first 
survey (e.g. Fukui et al. 1999; Mizuno et al. 2001; Fukui et al. 
2008). The observed grid spacing was 2\ corresponding to about 
30 and 35 pc at the distance of the LMC and SMC, while the 
half-power beam width was 2.6' at 115 GHz. 

We used the CO maps of the Magellanic Clouds to construct 
an integrated intensity map {Wco)^ integrating over the full visr 
range of the data (100 < visr < 400 kms"^for about 80 % of 
the data, while the remaining 20 % had a velocity range of about 
100 km s"^ covering the H i emitting regions.) 

2.2.4. data 

In order to estimate the free-free contribution to the millime- 
tre fluxes, we used the continuum- subtracted Hq, maps from the 
Southern H-Alpha Sky Survey Atlas (SHASSA, Gaustad et al. 
2001) centered on the Magellanic Clouds. 

2.2.5. FIR-Submm data 

We used the following FIR-submm ancillary data. 

- IRAS-IRIS (Improved Reprocessing of the IRAS Survey) 
100 |Ltm in order to constrain the dust temperature. The char- 
acteristics of this data, including the noise properties were 
taken from Miville-Deschenes & Lagache (2005) and are 
summarized in Table 1 

- WMAP 7yr data. The characteristics of this data, including 
the noise properties were taken from Bennett et al. (2003) 
and Jarosik et al. (2010) and are summarized in Table 1 



2.3. Additional Data processing 

2.3.1 . Common angular resolution and pixelisation 

For data already available in the HEALPix (Gorski et al. 2005) 
format (e.g., WMAP), we obtained the data from the Lambda web 
site (http://lambda.gsfc.nasa.gov/). For data not originally pre- 
sented in the HEALPix format, the ancillary data were brought 
to the HEALPix pixelisation, using a method where the surface 
of the intersection between each HEALPix and FITS pixel of 
the original data was used as a weight to regrid the data. The 
HEALPix resolution was chosen so as to match the Shannon 
sampling of the original data at resolution 6, with a HEALPix 
resolution set so that the pixel size is < 6/2 A. The ancillary 
data and the description of their processing will be presented 
in Paradis&et. al. (2011). 

All ancillary data were then smoothed to an appropriate res- 
olution by convolution with a Gaussian smoothing function with 
appropriate FWHM using the smoothing HEALPix function, 
and were brought to a pixel size matching the Shannon sampling 
of the final resolution. 



2.3.2. CMB subtraction 



Residual in LMC/SMC simulations at 100 GHz 




Patch size (degrees) 



Fig. 2. Error due to the ILC CMB subtraction derived from 
Monte-Carlo simulations for both LMC and SMC at 100 GHz, 
as a function of the patch size. The values shown are the diff'er- 
ence between the recovered and the input CMB divided by the 
simulated LMC and SMC, both integrated in a 4° ring. 



The upper panel in Fig. 3 shows the total intensity maps ob- 
served toward the LMC and SMC before CMB subtraction in 
the LFI 70.3 GHz chanel. The amplitude of CMB fluctuations 
is of the order of that of the diff'use emission of the galaxies and 
it is clear that, in order to study the integrated SED of the LMC 
and SMC galaxies, an efl&cient CMB subtraction has to be per- 
formed. 

The standard CMB-subtracted maps produced by the Data 
Processing Center (DPC) (The Planck Collaboration 2011g) 
were not used in this analysis. They were processed by a Needlet 
Internal Linear Combination (The Planck Collaboration 2011g) 
(NILC) that left a significant amount of foreground emission in 
the CMB estimate towards the LMC and the SMC. 

For this reason we have subtracted an estimate of the CMB 
optimized locally for the LMC and SMC regions, as described 
below. This CMB component was reconstructed through a clas- 
sical Internal Linear Combination (ILC) by means of Lagrange 
multipHers (Eriksen et al. 2004). 12° x 12° patches around the 
LMC and SMC were extracted from the HFI CMB frequency 
channels maps (100, 143, 217and 353 GHz) reduced to a com- 
mon resolution (10 arcmin), in units of Kcmb- The CMB com- 
ponent obtained on these patches clearly contains less LMC and 
SMC residual than the standard HFI DPC CMB and therefore 
will less aff'ect the SED determination. The lower panel in Fig. 3 
shows the maps after CMB subtraction in the LFI 70.3 GHz 
chanel. 

We performed Monte-Carlo simulations in order to estimate 
the error induced on the SED by our CMB removal. The LMC 
and SMC were simulated as a sum of two correlated compo- 
nents. The first component spatial template is the IRAS -IRIS 
100 |Ltm map. The second component spatial template is the 
545 GHz LMC or SMC Planck map. Their correlation coeffi- 
cients are 83% for the LMC and 92% for the SMC. We normal- 
ized their fluxes inside a 4° ring to the value of a typical LMC 
or SMC dust component and a typical millimetre excess at each 
HFI CMB frequency, respectively. 200 independent realizations 
of a WMAP 7yr best fit CMB (Komatsu et al. 2010) and of nomi- 
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Fig. 3. SMC (left) and LMC (right) total intensity maps before 
the 13.01 arcmin resolution. 



nal inhomogeneous HFI white noise were added to the synthetic 
LMC and SMC on patches of varying sizes. For each of these 
simulations, for each size, our ILC is performed and compared 
to the input CMB. We estimate the error due to the CMB sub- 
traction in both the LMC and the SMC by comparing the resid- 
ual in the CMB map to the emission of the two components at 
a given frequency. Results at 100 GHz are displayed in Fig. 2. 
For both the LMC and the SMC, the error increases at small and 
large patch sizes and an optimal patch size with respect to the 
CMB subtraction can be found at 5° for the LMC and 13° for 
the SMC. This behavior is due to the fact that the narrower the 
patch, the lower the contribution of the CMB to the total vari- 
ance which is to be minimized in the ILC. On the other hand, 
when their sizes increase, patches include foreground emission 
which is uncorrelated with the galaxies and has a different spec- 
trum, and thus both small and large patches contribute to the total 
variance which must be minimized. 

For the 12° x 12° patches used in the following analysis, we 
estimate the error due to CMB subtraction as 10.38 jliKcmb and 
28.2 illKcmb for the LMC and SMC respectively. In terms of the 
fraction of the total galaxy brightness, these correspond to 8.1, 
6.2, 2.3 and 0.3 % for the LMC and 12.1, 10.7, 6.4 and 1.4 % for 
the SMC at 100, 143, 217and 353 GHz, respectively. 




06*^00™ 05^^45'^ 30 15 00 04^^45™ 30 




«^ooo (hr) 



(top) and after (bottom) CMB subtraction in the 70.3 GHz band at 



2.3.3. Galactic foreground subtraction 

The Milky Way (MW) emission is non-negligible compared to 
the emission of the LMC and SMC. We remove this contribu- 
tion at all wavelengths using the MW H i template described in 
Sec. 2.2.2. We first computed the correlation between all FIR- 
submm data with this template, in the region where both the 
MW H I template and the CMB estimate are available, but ex- 
cluding a circular region centered on each galaxy (centre coor- 
dinates taken as 0^2000 = 05hl5m30s, ^2000 = -68°30'10" and 
0^2000 = 01h04ml6s, ^2000 = -72°5r36" for the LMC and 
SMC respectively) with radius 4.09° and 2.38° for the LMC 
and SMC respectively. The spectral distribution of this correla- 
tion factor, taken to represent the SED of the MW foreground is 
shown in Fig. 4. The SED is compared to that of the high galac- 
tic latitude reference region used in The Planck Collaboration 
(2011s) in Fig. 4. It can seen that the two SEDs are similar, al- 
though the MW foreground towards the LMC and SMC appears 
slightly colder and has a relatively stronger non-thermal compo- 
nent in the millimetre. We subtracted the MW foreground from 
all data using this SED multiplied by the MW H i template over 
the full map extent. Note that the median foreground Hi inte- 
grated intensities over the LMC and SMC are ^ 297Kkms"^ 
and ^ 172Kkms"\ which correspond to a brightnesses of 
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Fig. 4. Average foreground SED in the direction of the LMC and 
SMC (lower curve) including IRAS - IRIS (light blue), Planck- 
HFI (red), Planck-LFI (green) and WMAP data (dark blue), 
compared to the SED of the high latitude low column-density 
MW SED (upper curve) derived in The Planck Collaboration 
(201 Is), which has been scaled up by a factor of 10 for clarity. 



S.OMJysr"^ and l.TMJysr"^ respectively at 857 GHz. This is 
of the same order as the average brightness of the galaxies at 
that frequency (see Table 2). However, most of the MW emis- 
sion is canceled when subtracting a local background around the 
galaxies and the differential correction due to the spatial struc- 
ture of the MW foreground then accounts for about 0.4% and 
21% of the LMC and the SMC brightness respectively. 

3. Integrated SEDs 

The integrated SED of the LMC and SMC before and after CMB 
subtraction are shown in Fig. 5 (red and blue diamonds, respec- 
tively). They were computed by averaging values in the circular 
area around each galaxy defined in Sec. 2.3.3, and subtracting 
an estimate of the sky off' the source taken in an annulus with 
radius 1°. They are compared to the data taken from Israel et al. 
(2010) and Bot et al. (2010b), which were integrated over the 
same region. It can be seen that the ffux observed in the HFI and 
LFI bands is consistent with that obtained with previous studies 
in this wavelength range before subtraction of the CMB ffuctua- 
tions. A model of thermal dust, free-free and synchrotron emis- 
sion is fitted to the data as follows. As in Bot et al. (2010b), 
the thermal dust emission is adjusted according to the Draine & 
Li (2007) dust model. The free-free emission is deduced from 
the Hq, integrated ffux, using the expression from Hunt et al. 
(2004), assuming an electronic temperature Te = 10"^ K, the 
ratio of ionized helium to hydrogen nHe+/nJj = 0.087, and no 
extinction. The synchrotron emission was fftted to the radio data 
from the literature as in Israel et al. (2010). The combined dust, 
free-free and synchrotron emission is shown by the blue line. 
It can be seen from Fig. 5 that the SED after subtraction of the 
CMB ffuctuations is in fact compatible with no millimeter emis- 
sion excess for the LMC, for this particular model. However, 
the CMB -subtracted SED still shows a signiffcant excess for the 
SMC, with most data points above A ^ I mm being in excess 
over the model by more than 5cr, leading to an overall signiff- 
cance of the excess of about SOcr. 

The CMB subtraction removes part of the millimetre excess 
in both galaxies. This shows that CMB ffuctuations behind the 



Table 3. Dust temperature, JS values and fit reduced for vari- 
ous methods experimented to derive the temperature maps. The 
values listed are median values in the SED integration region for 
each galaxy. 



Method 


Td + Arc 


/3±A/3 






[K] 






LMC: 


freeyS 


21.0+1.9 


1.48+0.25 


1.91 


fixed 


20.7+1.7 


1.5t 


1.73 


fixed 


19.2+1.6 


1.8 


1.63 


fixed 


18.3+1.6 


2.0 


1.57 


dustem 


17.7+1.6 




1.60 


SMC: 


freeyg 


22.3+2.3 


1.21+0.27 


2.28 


fixed yg 


21.6+1.9 


1.2t 


1.90 


fixed yg 


18.8+1.7 


1.8 


9.94 


fixed yg 


17.9+1.6 


2.0 


14.0 


dustem 


17.3+1.6 




12.66 



fixed J3 model used in this paper. 



LMC and the SMC average out to a small but positive contri- 
bution when integrated over the extent of the galaxies. We note 
that an excess emission remains in the SMC independently of 
the dust model used and the assumptions made on the free-free 
or synchrotron emission. The shape and intensity of the mil- 
limetre excess can change accordingly, but we could not find a 
solution where the SMC SED is explained purely by thermal 
dust emission, free-free and synchrotron. We emphasize also 
that the dust model, used to reproduce the dust emission up to 
the sub-millimetre, assumes components heated by a radiation 
field 10 times lower than the solar neighborhood radiation field. 
Compared to other nearby galaxies, this result is rather extreme 
(Draine et al. 2007; Bot et al. 2010b). 

The integrated SED of the LMC and SMC after CMB and 
foreground subtraction are compared in Fig. 6. These SED were 
computed in the same integration region as those in Fig. 5. The 
comparison in Fig. 6 shows that the SMC SED is ff alter than the 
LMC one in the submm, while the SEDs when normalized in the 
FIR, reconcile above 10 mm, where the emission is presumably 
dominated by free-free and/or spinning dust (see Sec. 5). 

The SEDs at various stages of the background and fore- 
ground subtraction are given in Table 2 for the LMC and SMC. 
The uncertainties given include the contribution from the data 
variance combined for the integration region, the data variance 
combined for the background region, and the absolute calibra- 
tion uncertainties, using the values given in Table 1. They also 
include the noise resulting from the background (CMB) subtrac- 
tion and from the foreground MW subtraction and the free-free 
removal. All uncertainty contributions were added quadratically. 

4. Dust temperature and emissivity 

4.1. Temperature determination 

As shown by several previous studies (e.g. Reach et al. 
1995; Finkbeiner et al. 1999; Paradis et al. 2009; The Planck 
Collaboration 2011s,t), the dust emissivity spectrum in our 
Galaxy cannot be represented by a single dust emissivity in- 
dex JS over the full FIR-submm domain. The data available indi- 
cate that is usually stronger (steeper emissivity) in the FIR and 
lower (ffatter emissivity) in the submm, with a transition around 
500 |Ltm (see Paradis et al. 2009; The Planck Collaboration 
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Fig. 5. Integrated SEDs of the LMC (left) and SMC (right) before and after CMB subtraction. The black points and model are taken 
from Bot et al. (2010b). The red symbols show the SEDs derived from the DIRBE, IRAS, and WMAP data before CMB subtraction. 
The blue symbols show the same after CMB subtraction. 



Table 2. LMC (column 2-4) and SMC (column 5-7) SEDs averaged in a circular region for each galaxy. The integration region is 
centered on CK2000 = 05hl5m30s, ^2000 = -68°30'10", with radius Rlmc = 4.09° faor the LMC and centered on 0^2000 = 01h04ml6s, 
<^2000 = -72°5r36", with radius Rsmc = 2.38° for the SMC. A common background was subtracted in a 1° annulus around 
this region. Brightness values are in MJy sr"^ in the vly = cste flux convention. The last two lines give the average H^, emission 
(in Rayleigh), Galactic HI emission (in Kkms"^) in the same area. The table lists the total SED (f^^), the CMB subtracted SED 
^^noCMB^ and the CMB, MW foreground and free-free subtracted SED (7^^^) and their associated Icr uncertainties. 



A 


jtot 


jnoCMB 


jsub 


jtot 


jnoCMB 


jsub 


[|^m] 


[MJy sr-i] 


[MJy sr-1] 


[MJy sr-i] 


[MJy sr-i] 


[MJy sr-1] 


[MJysr-1] 


IRAS: 


12 


(2.36+0.12)xl0-^ 


(2.36+0.12)xl0-^ 


(2.34+0.13)xl0-^ 


(1.37+0.11)xlO-' 


(1.37+0.11)xlO-' 


(-2.79+0.27)xl0-' 


25 


(6.28+0.95)xl0-i 


(6.28+0.95)xlO-i 


(6.25+0.96)xl0-i 


(7.01 + 1. 10)xl0-2 


(7.01+1. 10)xl0-2 


(5.30+0.98)xlO-2 


60 


5.79+0.60 


5.79+0.60 


5.78+0.61 


1.52+0.16 


1.52+0.16 


1.56+0.18 


100 


(1.13+0.15)xlOi 


(1.13+0.15)xlOi 


(1.13+0.15)xlOi 


2.56+0.35 


2.56+0.35 


2.82+0.39 


Planck: 


349.82 


4.97+0.35 


4.97+0.35 


4.96+0.35 


1.18+0.08 


1.18+0.08 


1.46+0.11 


550.08 


1.79+0.13 


1.79+0.13 


1.78+0.13 


(5.02+0.37)xl0-i 


(4.99+0.37)xl0-i 


(5.96+0.45)xl0-i 


849.27 


(4.61+0.10)xl0-i 


(4.58+0.10)xlO-i 


(4.51+0.10)xl0-i 


(1.63+0.04)xl0-i 


(1.47+0.04)xl0-i 


(1.69+0.05)xl0-i 


1381.5 


(1.16+0.02)xl0-i 


(1.10+0.03)xlO-i 


(1.04+0.03)xl0-i 


(7.12+0.18)xl0-2 


(4.30+0.19)xl0-2 


(4.63+0.20)xl0-2 


2096.4 


(3.56+0.08)xlO-2 


(3.13+0.10)xlO-2 


(2.5 l+0.09)x 10-2 


(3.71+0.10)xl0-2 


(1.54+0.11)xl0-2 


(1.54+0.11)xl0-2 


2997.9 


(1.81+0.05)xlO-2 


(1.53+0.06)xlO-2 


(8.96+0.48)xlO-3 


(2.25+0.07)xl0-2 


(8.29+0.71)xl0-3 


(7.09+0.69)xl0-3 


4285.7 


(1.14+0.07)xl0-2 


(9.87+0.68)xlO-3 


(3.34+0.36)xlO-3 


(1.24+0.09)xl0-2 


(4.55+0.61)xlO-3 


(3.16+0.54)xl0-3 


6818.2 


(9.24+0.52)xl0-3 


(8.57+0.51)xlO-3 


(1.72+0.17)xl0-3 


(6.37+0.49)xl0-3 


(2.94+0.34)xl0-3 


(1.57+0.28)xl0-3 


10000 


(8.43+0.61)xl0-3 


(8.17+0.60)xlO-3 


(1.05+0.25)xlO-3 


(3.87+0.55)xlO-3 


(2.53+0.49)xl0-3 


(1.10+0.41)xl0-3 


WMAP: 


3200 


(1.61+0.04)xl0-^ 


(1.36+0.05)xlO-' 


(7.23+0.43)xl0--^ 


(2.10+0.08)xl0-' 


(8.28+0.80)xlO-^ 


(7.07+0.79)xl0-^ 


4900 


(1.04+0.02)xl0-2 


(9.13+0.23)xlO-3 


(2.51+0.17)xlO-3 


(1.06+0.03)xl0-2 


(4.32+0.33)xl0-3 


(2.92+0.31)xl0-3 


7300 


(9.15+0.13)xlO-3 


(8.56+0.15)xlO-3 


(1.66+0.09)xlO-3 


(6.12+0.15)xl0-3 


(3.09+0.16)xlO-3 


(1.67+0.15)xl0-3 


9100 


(9.38+0.12)xl0-3 


(8.98+0.13)xlO-3 


(1.93+0.06)xlO-3 


(5.10+0.11)xlO-3 


(3.04+0.11)xl0-3 


(1.63+0.10)xlO-3 


13000 


(9.38+0.11)xlO-3 


(9.19+0.12)xl0-3 


(1.88+0.05)xlO-3 


(3.62+0.08)xl0-3 


(2.52+0.09)xl0-3 


(1.07+0.07)xl0-3 


H« (R): 




(1.93+0.19)xl0^ 






7.00+0.72 








Kkms-^: 




1.01+0.13 






(-2.57+0.26) X 10+^ 







201 It). This is likely to be the case also for the LMC and the 
SMC. The dust temperature derived will depend on the assump- 
tions made about JS, since these two parameters are somewhat 
degenerate in space. Note that some of the past attempts at 
constraining the equilibrium temperature of the large grains in 
the LMC used the IRAS 60 jam emission. Emission at 60 \im 
is highly contaminated by out-of-equilibrium emission from 



Very Small Grans (VSGs) and this is even more the case in the 
Magellanic Clouds, due to the presence of the 70 |Ltm excess (Bot 
et al. 2004; Bernard et al. 2008). Combining the IRAS 60 |Ltm and 
100 |Ltm data therefore strongly over-estimates the temperature 
and accordingly under-estimates the abundances of all types of 
dust particles. A good sampling at frequencies dominated by Big 
Grain (BG) emission became possible using the combination of 
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Fig. 6. Integrated SEDs of the LMC (solid) and SMC (dashed) 
after CMB and galactic foreground subtraction, including data 
from Planck-HFI (red), Planck-LFI (green), IRAS - IRIS (light 
blue) and WMAP (dark blue). The uncertainties shown are ±3cr. 
The SMC SED was scaled by a factor 4 providing normalization 
at 100 |Ltm. 



the IRAS and Spitzer data (Leroy et al. 2007; Bolatto et al. 2007; 
Bernard et al. 2008; Sandstrom et al. 2010). 

As dust temperature is best derived from the FIR data, we 
limit the range of frequencies used in the determination to 
the FIR, which limits the impact of potential changes of the 
dust emissivity index jS with frequency. In the determination 
of the dust temperature (T^d), we used the IRAS -IRIS 100 jam 
map and the two highest HFI frequencies at 857 and 545 GHz. 
Temperature maps were derived at the common resolution of the 
3 bands used (5') and at lower resolution for further analysis. In 
each case, the emission was computed in the photometric chan- 
nels of the instruments used (IRAS, Planck, and WMAP), includ- 
ing the color corrections, using the actual transmission profiles 
for each instrument, and following the flux convention descrip- 
tion given in the respective explanatory supplements. 

In order to derive the thermal dust temperature, we use the 
same strategy as described in The Planck Collaboration (201 Is). 
To minimise computation time, the predictions of the model 
were tabulated for a large set of parameters (T^d, ft). For each 
map pixel, the was computed for each entry of the table and 
the shape of the distribution around the minimum value was 
used to derive the uncertainty on the free parameters. This in- 
cludes the efl'ect of the data variance (Tjj and the absolute uncer- 
tainties. 

We explored several options for deriving maps of the appar- 
ent dust temperature, which are summarized in Table 3. 

We first fitted each pixel of the maps with a modified black 
body of the form ly oc v^By(TD) in the above spectral range 
(method referred to as "free yS" in Table 3). The median values 
of Td and jS derived using this method are given in the first 
line of Table 3. This led to median values of y^LMC - 1-5 and 
y^SMC - 1.2 for the LMC and SMC respectively. Note that the 
value for the LMC is consistent with that derived using a combi- 
nation of the IRAS and Herschel data by Gordon et al. (2010). 
However, inspection of the corresponding maps shows corre- 
lated variations of the two parameters which are also correlated 
with the noise level in the maps. This suggests that the spuri- 
ous values probably originate from the correlation in parameter 



space and the presence of noise in the data, particularly in low 
brightness regions of the maps. 

We then performed fits of the FIR emission using the fixed 
P values, with the p values given above for the two galaxies 
(method referred to as "fixed and marked with a Mn Table 3). 
Although the median reduced was slightly higher than for the 
"free yS" method, the temperature maps showed fewer spurious 
values, in particular in low brightness regions. This resulted in 
a more coherent distribution of the temperature values than in 
the "free yS" case. Since we later used the temperature maps to 
investigate the spectral distribution of the dust optical depth, and 
the dust temperature is a source of uncertainty, we adopted the 
"fixed yS" method maps in what follows. The corresponding tem- 
perature maps for the LMC and SMC are shown in Fig. 7. 

The calculations were also carried out fox p = 1.8, which is 
the average Galactic value (The Planck Collaboration 2011s,t), 
in order to be able to derive dust emissivity values under the 
same assumption as in the MW. We also experimented using the 
"standard" p=2 value. 

4.2. Angular distribution of dust temperature 
4.2.1. LMC 

The temperature map derived here for the LMC shows a similar 
distribution to the one derived from IRAS and Spitzer data by 
Bernard et al. (2008). The highest temperatures are observed to- 
ward 30-Dor and are of the order of 24 K. The large scale distri- 
bution shows the existence of an inner warm arm, which follows 
the distribution of massive star formation as traced by known 
H II regions. Figure 8 shows that the warm dust at Td > 20 K in 
the LMC is reasonably well correlated with H^, emission, indi- 
cating that it is heated by the increased radiation field induced 
by massive star formation. Two regions departing from the cor- 
relation are visible to the SW of the 30-Dor at 0^2000 =05h40m, 
^2ooo=-70°30' and Qf20oo-05MOm, ^2000 =-72° 30'. The first re- 
gion was already identified in Bernard et al. (2008) who showed 
that this low column density region was unexpectedly warm 
given that no star formation is taking place in this area. They 
proposed that this could be an artifact of the IRAS -IRIS 100 |Ltm 
used in the analysis, possibly due to IRAS gain variations upon 
passage on the bright 30-Dor source. However, the analysis car- 
ried out here and the comparison between temperature maps ob- 
tained with and without MW background subtraction favors a 
biased temperature due to low surface brightness. This points to 
the necessity of correcting for temperature bias at low bright- 
ness. The second region is close to the edge of the map and is 
probably also aff'ected by low surface brightness. 

The Planck data however reveals cold regions in the outer 
regions of the LMC, which were not seen using the Spitzer data. 
This is essentially due to the large coverage of the Planck data 
with respect to the limited region, which was covered by the 
Spitzer data. The south of the LMC exhibits a string of cold 
regions with Td < 20 K. These regions correspond without ex- 
ception to known molecular clouds when they fall in the region 
covered by the CO survey. They also correspond to peaks of the 
dust optical depth as derived in the following sections. Similarly, 
a set of cold regions exist at the nortwest periphery of the LMC. 
The comparison with the distribution of CO clouds is shown in 
Fig. 8. As already noticed in Bernard et al. (2008), there is no 
systematic correlation between cold dust and the presence of 
molecular material, at least in the inner regions of the LMC. 
This was confirmed by the statistical characteristics of the IR 
properties of LMC molecular clouds established by Paradis et al. 
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Fig. 7. Upper panels: Dust temperature maps for the LMC (left) and SMC (right) computed from the foreground subtracted maps 
using the IRAS-IRIS 100 jam, HFI 857 and 545 GHz maps, using a fixed JSlmc = 1.5 and JSsmc = 1-2. Lower panels: relative 
uncertainties on the dust temperature at the same resolution, expressed as percentages. 



(2010), which showed no systematic trend for molecular regions 
to be colder than their surrounding neutral material. However, 
toward the outer regions of the LMC, molecular clouds appear 
systematically to show a decrease in the dust temperature. This 
diff'erence may be due to the absence of star formation activity 
in the outer regions and/or to less mixing along the line of sight. 



Gordon et al. (2010) derived a dust temperature map from 
the Herschel data (HERITAGE program) for a fraction of the 
LMC observed during the Science Demonstration Phase (SDP). 
They used the 100 |Ltm to 350 jam IRAS and Herschel bands to 
constrain the temperature. They found that the temperatures de- 
rived this way only diff'er from those derived by Bernard et al. 
(2008) by up to 10%. The SDP field covered an elongated re- 
gion across the LMC roughly oriented north-south. This strip 
crossed the warm inner arm, which is also clearly seen in their 
map. However, we stress that, in their study, a gradient was re- 
moved along the strip based on the values at the edges. The fact 
that we detect significantly colder than average dust in the outer 
cold arm of the LMC underlines the need to take those varia- 
tions into account when subtracting background emission in the 
Herschel data. 



4.2.2. SMC 



The temperature map derived for the SMC shows more 
moderate temperature variations than those observed in the 
LMC but peaks corresponding to well known HII regions 
are clearly identified, and the overall dust temperature cor- 
relates with Hq, emission, as shown in Fig. 9. In particu- 
lar, the massive star forming region SMC-N66 (Henize 1956) 
(a^2000=00h59m27.40s, ^20oo=-72°10ai") corresponds to the 
highest temperature (~ 20 K) in the SMC. Other well known star 
forming regions like SMC-N83/84 (01hl4m21.0s, -73° 1742"), 
N81 (01h09ml3.6s, -73° 11^41"), N88/89/90 (in the wing: 
01h24m08.1s, -73°08'55") and DEM S54 (at the center of the 
main bar: 00h50m25.9s, -72°53'10") appear as temperature 
peaks compared to the surroundings. In contrast, the infrared- 
bright star forming region N76 (located at the northeast of N66: 
01h04m01.2s, -72°0r52") and southwest star-forming com- 
plex do not stand out. This could be due to the presence of an 
extended warm component that we observe in the main bar and 
that seems spatially related to the diff'use Hq, emission. 
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Fig. 8. Comparison between the dust temperature map of the 
LMC with Hq, (top) and CO emission (bottom). The CO contours 
are at 0.5, 2, 4 and 10 Kkms"^ contours are at 1, 10, 50, 
100, 500 and 1000 Rayleigh. The thick line shows the edge of 
the available CO surveys. 



Fig. 9. Comparison between the dust temperature map of the 
SMC with Hq, (top) and CO emission (bottom). The CO contours 
are at 0.5, 1 and 1.5 Kkms"^ Hq, contours are at 1, 5, 30 and 
100 Rayleigh. The thick line shows the edge of the available CO 
surveys. 



4.3. Optical depth determination 
Optical depth is derived using: 

T(V) : 



(1) 



where By is the Planck function. We used resolution matched 
maps of Td and ly and derived r maps at the various resolutions 
of the data used here. The uncertainty on r (At) is computed 
accordingly as: 



AT(y) = T 



n ^\ ST BAT^) 



2^1/2 



(2) 



The optical depth and optical depth uncertainty maps derived 
at 857 GHz are shown in Fig. 10 and Fig. 1 1 for the LMC and the 
SMC respectively. 



5. Discussion 

5.1. The millimetre excess 

Bot et al. (2010b) (see Fig. 5) found that the integrated SED of 
the LMC and more noticeably of the SMC showed excess mil- 
limetre emission with respect to a dust and free-free/synchrotron 
model based respectively on the FIR and the radio data avail- 
able. They investigated several causes for this excess. They no- 
ticed that the excess had precisely the colors of CMB fluctu- 
ations. They performed simulations by placing the LMC and 
SMC structure at various positions over a CMB simulated sky, 
and concluded that a CMB origin was unlikely, but not excluded. 
Other proposed origins for the excess included the presence of 
very cold dust, spinning dust or large modifications of the optical 
properties of thermal dust. 
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Fig. 10. Upper panel: Map of the dust optical depths of the LMC 
diiHFI 217 GHz. Units are Wxt. Lower panel: Map of the dust 
optical depth relative uncertainty of the LMC at HFI 217 GHz 
in percent. Black pixels in the maps are masked and have relative 
uncertainties larger than 50%. 




Fig. 11. Upper panel: Map of the dust optical depths of the SMC 
diiHFI 217 GHz. Units are Wxr. Lower panel: Map of the dust 
optical depth relative uncertainty at HFI 217 GHz in percent. 
Black pixels in the maps are masked and have relative uncer- 
tainties larger than 50%. 



Very cold dust (Td - 5 - 7 K) has been advocated to explain 
the flattening of the millimetre emission observed in more dis- 
tant low metallicity galaxies (e.g. Galliano et al. 2005). However, 
this usually led to very large masses, and the existence of such 
very cold dust remains controversial and difficult to understand 
in low metallicity systems where the stronger star formation rate 
and the lower dust abundances prevents eflftcient screening from 
UV photons. 

Bot et al. (2010b) applied spinning dust models to fit the 
SMC and LMC SEDs and found a plausible match. However, 
since the observed excess peaked at a significantly higher fre- 
quency than observed for spinning dust in other regions (e.g. The 
Planck Collaboration 201 Iv), their fit required extreme density 
and excitation conditions for the small dust particles. 

Using the Two Level System model by Meny et al. (2007) for 
the long wavelength emission of amorphous solids also proved 



plausible for the LMC, but a convincing fit could not be found 
for the SMC, essentially because the model could not reproduce 
the shape of the excess. 

The study carried out here, which takes advantage of the 
Planck measurements to constrain the CMB foreground fluctu- 
ations towards the two galaxies, shows that part of the excess 
observed toward the SMC cannot be accounted for by the fluc- 
tuations of the background CMB, as discussed in Sec. 2.3.2, but 
the intensity of the excess has been greatly reduced compared to 
the one found by Bot et al. (2010b). 

To give more insight into this excess, we built a map to trace 
the spatial distribution of the excess emission in the SMC. To do 
so, we applied the fitting procedure performed for the integrated 
SEDs of the Magellanic Cloud, to each point of the SMC at the 
angular resolution of the LFI smallest frequency channel. The 
SEDs were build at each point using the CMB and foreground- 
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Fig. 12. Map of the millimetre excess in the SMC at 3 mm com- 
puted at the resolution of the LFI -28.5 GHz channel, after free- 
free subtraction. The contours show the Hq, distribution. 



Gordon et al. (2010) found hints of such excess emission 
around 500 jam in the LMC Herschel data, but this remained 
within the current calibration uncertainties of the Herschel 
Spectral and Photometric Imaging Receiver (SPIRE ) instru- 
ment. We also observe that the 500 jam optical depth of the LMC 
is slightly higher (by 8.4%) than the power law extrapolation 
shown in Fig. 13, but this is only marginally larger than the Icr 
uncertainty at this wavelength. In any case, data points at longer 
wavelengths do not support the existence of excess emission in 
addition to a/5 = 1.5 power law for the LMC. 



subtracted data. To model the SEDs, we assumed that, at each 
point, the SEDs are dominated by free-free and dust emission up 
to the longest wavelengths covered by Planck, following what 
is observed for the integrated emission of the SMC (see Fig. 5). 
The free-free component at each wavelength was extrapolated 
from the Hq, emission (as in Sect. 3), assuming that the extinc- 
tion in the SMC is negligible. The thermal dust emission was 
then fitted to the data points at /Iref < 550 |Ltm, using the Draine 
et al. (2007) dust model. The millimeter excess was then de- 
fined as the diff'erence between the data and the dust-and-free- 
free model. The resulting spatial distribution of the excess emis- 
sion at 100 GHz (3 mm) is shown in Fig. 12. This shows that the 
peak of the excess is located at the southwest tip of the bar. This 
region also corresponds to the maximum of the optical depth de- 
rived from the FIR and shown in Fig. 1 1 and the overall excess 
spatial distribution is consistent with being proportional to the 
dust column density. 

5.2. FIR dust emissivity 

The wavelength dependence of the average dust optical depth in 
the LMC and SMC is shown in Fig. 13. The spectral index of the 
emissivity in the FIR is consistent with j3 = 1.5 and j3 = 1.2 for 
the LMC and the SMC respectively. This value for the LMC 
is consistent with findings by Gordon et al. (2010) using the 
Herschel data. We see no hint for for a change of the spectral 
index with wavelength for the LMC. In contrast, the SMC SED 
clearly flattens at A > 800 |Ltm to reach extremely flat j3 values 
around A = 3 mm. The dust emissivities as interpolated using 
the power laws shown in Fig. 13 can be compared to the ref- 
erence value for the solar neighborhood of t/Nh = 10"^^ cm^ 
at 250|Ltm (Boulanger et al. 1996). This comparison indicates 
lower than solar dust abundances for the two galaxies, by about 
1/2.4 and 1/13 solar respectively. This is in rough agreement with 
metallicity for the LMC and significantly lower than metallicity 
for the SMC. 



Fig. 13. Average dust optical depth of the LMC (black) and 
SMC (blue) obtained from the CMB, foreground and free-free 
subtracted SEDs. The average is taken in the same regions as for 
Fig. 5. The SMC SED was normalized (multiplied by 6.31) to 
that of the LMC at 100 |Ltm. The square symbol (green) shows the 
reference value for the solar neighborhood by Boulanger et al. 
(1996). The dashed lines show r oc y^-^ and r oc y^-^ normalized 
at 100 |Ltm. Error bars are ±3-cr. 



5.3. Possible interpretation of the excess 

In this section, we attempt to fit the SED of the SMC using vari- 
ous models which have been proposed to explain excess submm 
emission in addition to a single dust emissivity power-law for 
thermal dust. Figure 14 shows several fits to the FIR-Submm 
SED of the SMC. 

The first fit uses the Finkbeiner et al. (1999) model p (see 
their Table 3), which was designed to explain the flatter than ex- 
pected emission spectrum of our Galaxy as observed by FIRAS. 
Here we use this model to assess the possibility that the mil- 
limeter excess is due to very cold dust. The model was fit- 
ted using the same JS values as in Finkbeiner et al. (1999) for 
the two dust components O^warm = 2.6, y^coid = 1-5). It also 
assumes the same type of relationship between the cold and 
warm dust temperatures, which in the model reflects the fact 
that both dust species are subjected to the same radiation field. 
However, we allow the IR/optical opacity ratio (qcoid/qwarm) to 
vary, so that the ratio between the warm and cold dust temper- 
atures is free to vary. We also leave the cold dust component 
abundance (fcoid) and the warm dust temperature (Twarm) as free 
parameters. The best fit values are found for Twarm = 16.4 K 
qcoid/qwarm = 168.2, fcoid = 8.7 10"^. For thcsc parameters, the 
temperature of the cold component is Tcoid = 5.9 K. This cold 
temperature is needed to reproduce part of the millimetre excess 
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Fig. 14. Fit of the SMC SED using the Finkbeiner model (upper left), the TLS model (upper right), the spinning dust model (lower 
left) and a combination of TLS and spinning dust models (lower right). In the lower panels, the squares represent the flux in each 
band predicted by the best model (the blue line). 



but is much colder than that obtained by Finkbeiner et al. (1999) 
for the Galaxy (Tcoid = 9.6 K). This is obtained at the expense of 
strongly increasing the IR/optical opacity ratio for the cold com- 
ponent (qcoid/qwarm increased by a factor 15), which controls the 
Tcoid/Twarm ratio. The mass fraction of the cold component de- 
rived here, fcoid, is about 4 times lower than that derived for the 
MW in Finkbeiner et al. (1999), which compensates for the in- 
creased qcoid/qwarm- H is apparent that, despite invoking a much 
larger IR/optical opacity ratio for the cold particles, such a model 
has difficulties producing the submm excess above about X = 2 
mm. 

The second fit employs the Two-Level-System (TLS) model 
developed by Meny et al. (2007). The fit was obtained by mini- 
mizing;^^ in the range 100 jam < A < 5 mm against the follow- 
ing 3 parameters: Td the dust temperature, /cor the correlation 
length of defects in the material and A the density of TLS sites in 
the material composing the grains. The values derived for these 
parameters are 7^0= 18.9 K, /cor= 12.85 nm and A=7.678. The 
reduced for these parameters is =2.56. Compared to the 
best values found from Paradis (2007) for the MW (7^0=17.9 K, 
/cor= 12.85 nm, A=2.42), this indicates dust material with more 
TLS sites (more mechanical defects) but a similar defect corre- 
lation length compared with the dust dominating the MW emis- 
sion. Note that a grey-body fit of the same SED over the same 



frequency range leads to 7^0= 18.9 K and /3=0.14 and a reduced 
;^^=6.12, showing that the TLS model reproduces the data better 
than a single grey-body fit. 

The third fit uses spinning dust. In that case, the thermal 
dust emission was reproduced using the Draine & Li (2007) dust 
model. The remaining excess is then fitted with the spinning dust 
model described in Silsbee et al. (2010). Spinning dust emission 
has been proved to be sensitive to the neutral and ionized gas 
densities and to the size distribution of dust grains. The radi- 
ation field and the size distributions are taken to be the same 
as for the Draine & Li (2007) model. The relevant gas parame- 
ters are computed with CLOUDY (Ferland et al. 1998): we take 
the parameters from the optically thin zone of isochoric simula- 
tions. In both cases, we take into account the lower metallicity 
and dust grains abundance when compared to the MW. The elec- 
tric dipole moment distribution of grains is taken as in Draine 
& Lazarian (1998b)^. This fit is obtained by adding two com- 
ponents according to the PDR fraction inferred from the ther- 
mal emission model: a diffiise medium with ^H=30cm"^ and 
100% of the PAH mass expected from the IR modeling, and a 



^ Ysard & Verstraete (2010b) showed that it is in good agreement 
with the anomalous emission extracted from the WMAP data. 
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denser medium with ^H=5000cm"^ and 82% of the PAH mass 
expected. 

The fourth one is a combination of the TLS model and the 
spinning dust model presented above. 

It is clear from Fig. 14 that the model combining TLS and 
spinning dust is the only one of the proposed models giving a 
satisfactory fit to the millimetre excess observed in the SMC over 
the whole spectral range. In addition, it alleviates the need for us- 
ing more PAH than allowed by the NIR emission to produce the 
required level of spinning dust emission observed. Alternative 
models, however, cannot be excluded. For instance, very large 
grains (a > 30 jam) that would be efficient radiators (Rowan- 
Robinson 1992) could create emission in the millimetre wave- 
lengths without too much mass. Exploring this possibility would 
require specific modeling. This is beyond the scope of the cur- 
rent article but should be explored in future studies. 

6. Conclusions 

We assessed the existence and investigated the origin of mil- 
limetre excess emission in the LMC and the SMC using the 
Planck data. The integrated SED of the two galaxies before sub- 
traction of the foreground (Milky Way) and background (CMB 
fluctuations) emission are in good agreement with previous de- 
terminations. 

The background CMB contribution was first subtracted us- 
ing an Internal Linear Combination (ILC) method performed lo- 
cally around the two galaxies. The uncertainty of this contribu- 
tion was measured through a detailed Monte-Carlo simulation. 
We subtracted the foreground emission from the Milky Way us- 
ing a Galactic H i template and the proper dust emissivity derived 
in a region surrounding the two galaxies and dominated by MW 
emission. We also subtracted the free-free contribution from ion- 
ized gas in the galaxies, using the Hq, emission, taking advantage 
of the low extinction in those galaxies. The remaining emission 
of both galaxies correlates with the gas emission of the LMC and 
SMC. 

We showed that the excess previously reported in the LMC 
can be fully explained by CMB fluctuations. For the SMC, sub- 
tracting the CMB fluctuations decreases the intensity of the ex- 
cess but a significant millimetre emission above the expected 
thermal dust, free-free and synchrotron emission remains. 

We combined the Planck and IRAS-IRIS data at 100 |Ltm to 
produce thermal dust temperature and optical depth maps of the 
two galaxies. The LMC temperature map shows the presence of 
a warm inner arm already found with the Spitzer data, but also 
shows the existence of a previously unidentified cold outer arm. 
Several cold regions were found along this arm, some of which 
are associated to known molecular clouds. 

We used the dust optical depth maps to constrain the thermal 
dust emissivity spectral index {J3). The average spectral index in 
the FIR (/I < 500 jam) is found to be consistent with p = 1.5 
and p = 12 for the LMC and the SMC respectively. This is 
significantly flatter than what is observed in the Milky Way. The 
absolute values of the emissivities in the FIR, when compared 
to that in our solar neighborhood, are compatible with D/G mass 
ratio of 1/2.4 and 1/13 for the LMC and SMC. This is compatible 
with the metallicity for the LMC but significantly lower than 
metallicity for the SMC. In the submm, the LMC SED remains 
consistent with p - 1.5, while the SED of the SMC flattens even 
more. 

The spatial distribution of the mm excess in the SMC appears 
to follow the general pattern of the gas distribution. It therefore 
appears unlikely that the excess could originate from very cold 



dust. Indeed, this is confirmed by attempts to fit the SMC emis- 
sion SED with 2 dust component models, which led to poor fits. 
Alternative models, such as emission excess due to the amor- 
phous nature of large grains are likely to provide a natural ex- 
planation to the observed SED, although spinning dust is needed 
to explain the SED above /I =3 mm. 

Acknowledgements. A description of the Planck Collaboration and a list 
of its members can be found at http://www.rssd.esa.int/index.php? 
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